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Abstract 

Time dependence for barrier penetration is considered in the phase space. An 
asymptotic phase-space propagator for nonrelativistic scattering on a one - dimen- 
sional barrier is constructed. The propagator has a form universal for various initial 
state preparations and local potential barriers. It is manifestly causal and includes 
time-lag effects and quantum spreading. Specific features of quantum dynamics 
which disappear in the standard semi-classical approximation are revealed. The 
propagator may be applied to calculation of the final momentum and coordinate 
distributions, for particles transmitted through or reflected from the potential bar- 
rier, as well as for elucidating the tunneling time problem. 

PACS number. 03.65.Nk 

1 Introduction 

Observable properties of quantal systems, like energy levels and transition probabilities are 
mostly related to stationary states. Meanwhile, the time dependence of physical processes 
is also described by quantum theory and may be of considerable interest. An important 
class of effects is various barrier penetration (tunneling) processes. The transition prob- 
abilities are usually obtained by means of the time-independent (energy) methods, in 
particular, in the semi-classical approximation (see e. g. in Ref. 0]). As soon as one 
gets the complete solution in the energy representation, the time evolution is obtained 
straightforwardly, in principle, in terms of the inverse Laplace transform. However, the 
evaluation of the large-time asymptotics may be an intricate job. A source of the trouble 
is in the very statement of the problem. If we insist that the particle was on one side 
of the barrier in the beginning, the state cannot be described by the plane wave, which 
is the eigen-state of the momentum operator. Thus the initial energy is never free of an 
uncertainty. The uncertainty may be made smaller if the particle is de-localized in space, 
so one has to start far enough from the barrier, and to detect the result long enough after 
the start in order to be sure that the particle has left the potential domain completely. It 
is clear, however, that the problem needs a special theoretical analysis. 
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The time dependence of tunneling processes has been attracting attention for decades. 
A controversial question is that of the tunneling time and the effect of causality on the 
particle propagation New experimental techniques enable detailed measurements 

performed on a photon wave packet transmitted through an optical analog of a quan- 
tum potential barrier [ 17 1. There are some theoretical concerns on the validity of the 
semi-classical approximation in processes, like tunneling, which have no classical counter- 
parts, because the effect of quantal fluctuations must get a proper account. Therefore a 
consistent time-dependent formalism for tunneling is hardly redundant. 

Barrier penetration processes and their time dependence were investigated by numeri- 
cal, experimental, and analytic methods in a number of works, e.g.|18]-[^. The literature 
contains specific examples of barriers and wave-packet shapes. Our purpose was to con- 
sider a general case, with no assumptions on the initial state preparation, neither on the 
form of the (local) potential barrier. Barrier penetration is described sometimes by means 
of the imaginary time method p9|-|3^. In the present approach, the process is described 
in real space-time, even though we exploit analytical properties in the complex energy 
plane, especially for causality arguments. 

The method applied here is an investigation of the time evolution of the Wigner phase- 
space distribution [0]. In this way, we can consider any initial state, not necessarily pure, 
which is important for applications to experiments. Besides, cumbersome oscillations of 
the wave functions are not involved. The Wigner function was used successfully in many 
problems of quantum theory, and its properties where considered, e.g. in |^-@]. Car- 
ruthers and Zachariasen, in their review of quantum collision theory with phase-space 
distributions considered cross sections for scattering processes in three dimensions, 
inclusive reactions within a second quantization approach, inclusive multi-particle pro- 
duction processes in the ultrarelativistic domain and many other processes but discussed 
neither quantum jumps, nor barrier penetration. Various other approaches have been 
tried recently [H5|-E^ to describe quantum dynamics in terms of the Wigner function. In 



scattering problems, the time evolution of the Wigner function was considered mainly 
within the semi-classical approximation p0|-||53|| . As was shown by Berry |5^, the semi- 
classical approximation for stationary Wigner functions describing bound states is given 
by the Airy function, its spread from the classical 5-function being the first-order quantum 
effect. Propagation of wave packets was discussed in a number of works; it was shown 



in particular that the spreading in the coordinate is not a specific quantal effect ||55 



Within the semi-classical approximation the quantal features of the long-time evolution 
was attributed to the interference between amplitudes corresponding to different classi- 
cal paths[Q. Our point is to emphasize the difference between the genuine quantum 
dynamics and the semi-classical approximation for classically forbidden processes where 
no classical paths exist, so quantum dynamics is not reduced just to a smearing around 
classical paths, or to an interference between them. 

The Wigner function was applied to tunneling by Balazs and VorospTf for parabolic 
potential barrier. The equivalence of Wigner's integro-differential equation to Liouville's 
classical equation is in an apparent conflict to tunneling for that potential. The puz- 
zle is solved as soon as it is realized that the initial Wigner function cannot be chosen 
arbitrarily, if the potential does not vanish asymptotically. Physically available states 
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correspond to Wigner distributions extended in the momentum, so real classical trajec- 
tories transport the particle above the barrier. The Wigner function for the parabolic 
potential was constructed explicitly, and the transmission and reflection coefficients were 
obtained, within a qualitative picture of tunneling where the difference between classical 
and quantum mechanics lies mainly in the initial state preparation. 

For sufficiently broad incident momentum distributions, the classical paths enable the 
particle to overcome the barrier, and the semi-classical approach can be used indeed ||58|, 
pU[ . If the barrier potential is localized, the incident wave packet (or the Wigner function) 
can be prepared with an arbitrarily definite energy, and no classical paths would be 
responsible for the tunneling. In that situation the barrier penetration must be a result 
of an essential difference between the quantum and classical dynamics. 

Purely quantal effects and their role in scattering processes have been also considered in 
the phase space formalism. In the "Wigner-trajectory" approach each phase-space point in 
the initial Wigner distribution is propagating along a definite trajectory]^, ^ ^ . If the 
third and higher-orders derivatives of the potential vanish, the Wigner trajectories coincide 
with the classical paths. Otherwise, the Wigner trajectories are defined with a modified 
"quantum" potential and are not classical. For a system in an energy eigen-state (i.e. 
in the stationary barrier problem), the time-shift invariance implies that the trajectories 
are the "equi- Wigner curves" which are lines of constant values of the Wigner function. 
That approach seems rather problematic. The effective 'potential' may be singular, the 
Liouville theorem is violated, and quantum jumps [l62|, |63[| can hardly be included. In 
another approach, quantum corrections to classical dynamics are interpreted as finite 
momentum jumps between classical paths[^. Negative quasi-probabilities appear in the 



the calculation, which distinguish the quantum treatment from the classical theory. The 
phase-space points are smeared to finite domains and do not propagate along continuous 
trajectories. 



We consider the phase-space evolution kernel 1 64], which is the fundamental solution of 



the dynamical equation for the Wigner function. In classical theory, the evolution kernel 
is the fundamental solution of the Liouville equation and equals the 5-function restricted 
to classical trajectories. In the semi-classical approximation one can get the Airy function. 
For the barrier penetration, an explicit expression is obtained and it is shown that it is 
not reduced to the semi-classical approximation. From this point of view, the barrier 
penetration is an essentially quantal process. 

In Section 2 the phase-space propagator (the evolution kernel) is defined, and some 
of its properties are given. Section 3 shows the relation between the time evolution and 
the S'-matrix formalism in the momentum representation. The large-time asymptotics for 
the space-time propagator is derived in Section 4. The result is an integral representation 
in terms of scattering amplitudes. The exact result for the narrow potential barrier is 
presented in Section 5. The semi-classical approximation is considered in Section 6, and 
its vahdity is discussed. Some technical aspects are considered in Appendix: the accuracy 
of the large-time asymptotics, the transition probability for Gaussian states, and the exact 
result for the cosh~^ potential barrier. 

The units used in the paper are h = 1, and 2m = 1 for the particle mass. 
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2 The phase-space propagator 

In general, any quantum state is described by a density matrix p, which can be represented 
by its matrix elements, say, in the coordinate representation, {q'\p\q)-, or by its Weyl 
symbol (the Wigner function p4|) p(x), where x = 



piq,p) = / {q + i:\p\q -i)e~''"^d7]. 



For a given Hamiltonian the evolution operator Uif) = exp{—iHt) determines the 
time evolution of the state: pt = U{t)poU^{t). Respectively, the time evolution of the 
Wigner function is given by an integral kernel, i.e. the phase-space propagator, 

Ptix) = J Lt{x,xo)po{xo)dxo. (2) 

Here dx = dqdp, and Lf is a real dimensionless function which satisfies the following 
identities, 

Lt{q,p;qo,po) = L^tiqo, -Po] q, -p), 
/d.L,(.,.„) = l=/d.oi,(.,.,), 

L„,(.,.„)=/d.'L,(.,.')i,.(.',.o) (3) 

In order to handle experimental data, one may introduce an acceptance function ({x), 
characterizing the apparatus, so that the probability to detect the system described by 
p{x) would be 

wc^ = j ({x)p{x)dx. (4) 



Thus the probability to detect the system at a time t after it was prepared in the state 
given by po(^) is 

w^t) = / / C{x)Lt{x,xo)po{xo)dxodx. (5) 



The quasi-distributions p{x) and ({x) are normalizable and satisfy a set of conditions 
owing to the positive definiteness of the density matrix, for instance, 



[p{x)]^dx < 



p{x)dx 
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(6) 



Other conditions are more intricate. Qualitatively, the distributions cannot be localized 
to domains of areas less than 27rh by the order of magnitude. The functions are not 
necessarily positive everywhere, but the domains of negativity must be small enough. 

The phase-space propagator Lf is expressed in terms of the matrix elements of tJ(t), 
e.g. the coordinate (or momentum) propagators, as follows from (1), 

Lt{q,p; go,Po) = r T d^oe'^^^^^q - ||?7|go - + ||?7|go + ^) 

ZTT J-oo J-oo Z II A 

^ i XI^^ I'^J^^^'^'^'^V + \\U{t)\p, + y)(p - ||f>(t)bo - y). (7) 
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For any Hamiltonian, quadratic in x, pt{x) is the solution of the classical Liouville 
equation. In that case the classical equations of motion are linear, as well as the Heisenberg 
equations. The solution is linear: x — RtXo, where Rt is an x-independent matrix, and 
Lt{x;xo) — S{x — RtXo). In particular, for the nonrelativistic free motion one has 

Lt{q,p;qo,Po) ^ 5(p - po)5{q - vt - qo), (8) 

where v = p/m is the particle velocity. 

For non-linear systems, effects specific for quantum dynamics result in deviations of 
Lt from the 5-f unction. 



3 Time dependence in the momentum representation 

As soon as the Hamiltonian H has no explicit time dependence, the evolution operator 
can be written as the Laplace transform of the resolvent: 

U{t) = e-'*^ = — / Gee~'^'de , G, = {H -e)'^ . (9) 



2711 JVo 

where is the usual integration contour in the complex e-plane, running above the real 
axis. 

We shall consider a nonrelativistic particle of mass m = 1/2 scattered from a localized 
one-dimcnsional and time-independent potential barrier. The Hamiltonian is H = Hq + V ^ 
where Hq = p^ is the kinetic energy operator. The basis of the normalized momentum 
eigen-states \k) will be used, so 

Ho\k)^k^\k), {k\ko) = d{k-ko). (10) 
Introducing the transition operator T^, one has 

4 = - Gf tGf , G'f = (^o (11) 
The momentum propagator is now given by 

{k\U{t)\ko) ^ ^J'^ {k\G,\ko)e-^''de (12) 

1 f ^^^-ite {k\T^\ko) 



- S{k - ko)e-'''' - / 



{k^ — i'j — s){kQ — i'j — e) 

The infinitesimal positive quantity 7 is introduced here to specify the integral near the 
poles due to the free propagators. The integration contour Foo is deformed to run around 
the positive real the exponential vanishes in the lower half-plane and {k\T^\kQ) 

has singularities only for the real values of s corresponding to physical energy values. The 
kinematic poles are isolated, leaving an integral along the real positive axis 



{k\Uit)\ko)^eM-^t{e + K)] 



5{k-ko)-^{k\T,\ko) |e=fe2 (13) 



+ ^{k\T,\ko) |,=feg -Jtik.ko) 
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where ^ = -^{kl — k"^) and 

and ^0 is a threshold energy value. (With no bound states eg = 0, otherwise it is the 
lowest bound-state energy.) At two zeroes of the integrand denominator the integral is 
taken in the sense of its principal value. 

In order to calculate explicitly the time-dependent propagator for a given potential, 
the transition operator matrix elements for this potential should be known on and off the 
energy shell. For scattering problems, one needs the large-time asymptotics, where the 
matrix elements are reduced to the energy shell because of the known fact of the theory 
of distributions, 

lim e*«Ve = n^5{i). (15) 

t — ^oo 

The integral in (14) vanishes, as t ^ oo, because Im [(/clT^I/co)] is smooth and the integral 
is converging. In general, Jt = 0(t~^/^). (The proof is given in Appendix A.) 
The result is expressed in terms of a unitary 2x2 matrix S 

{k\U{t)\kQ) X e-'^''2K6{k^ - kl)S,,,, SS^ = /, (16) 

where k = \k\ = \kQ\, I is the unit 2x2 matrix and v = k/n = ±1. The S-matrix elements 
are related to elements of the T^-operator on the energy shell = e, 

nx ^ 

Suuoii^) = Suuo (z/kITeIz/qk). (17) 

These are the probability amplitudes for transmission through and reflection from the 
potential region. The amplitudes can be expressed in terms of two analytical functions 
a{K) and 6(k), 

^('^) = ^(_\ \ y \a\'-\b\' = l. (18) 

These functions are defined for Re k > by the asymptotics of the solution of the 
stationary Schrodinger equation, 

p'^y + V{q)y = K^y, 
y-iQ) ^ I -ita I'l ina ^ ^ 7"^' (19) 
The functions a{K,) and 6(k) have the analytical continuation to the left half-plane by 



a{-R) = a{K), b{-R) = h{K). (20) 

The analytical properties of these functions in the complex K-plane have been investigated 
previously |^5[. It was shown, in particular, that for any finite-range and positive potential 
they can be expressed in terms of two entire functions a{e) and /9(e), 

a(ft) = 1 - a(£)/2m, 6(k) = (3{e)/2m. (21) 
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Thus the only singularities of S{k) are poles due to zeroes of a{K.) which are all in the 
lower half of the K-plane. Moreover, a and (3 are bounded for Im k > 0, so that lim a — 1 
and lim 6 = 0, as |«;| — > oo in the upper half- plane. If the potential has an exponential 
decrease as q —>■ ±oo, the functions a and P may have infinite series of poles. Besides, 
for symmetric barriers, where V{—q) — V{q), one has a real P{e) and purely imaginary 

b{K) = -b{-K). 



4 The large-time asymptotics 

In order to get the large time asymptotics for the phase-space propagator, we set the 
amphtudes from (16) into Eq. (6). Using the following property of the 5-function 

2\k\6{k^ - kl) = S{k - ko) + S{k + ko), (22) 

one gets the following result 

LM,P]%,Po) 5{.p-Po)'T{pQ,rj^) + 6{p + po)Tl{po,r_) (23) 



2^ 
+-Re 

TT 



^(PO ^2i(qop-qpo)+Mppot 



a{po - p)a{po + p) 



Here r± = qQ + 2pot^q; i.e the differences between the free classical trajectory and actual 
positions of the transmitted particle (r+) and reflected from q = particle (r_). The first 
two terms describe transmission and reflection from the potential barrier. The third term 
represents an interference between transmitted and reflected waves and is responsible for 
quantum fluctuations at the barrier region. It is irrelevant for large t, if the wave packet 
was prepared in free space with a narrow momentum distribution. 

The functions representing the transmission and reflection probabilities are given by 
the Fourier integrals 

T{po,r+) = — / 24 

In J-oo a[^a + pQ)a{^a - poj 

27T J-oo al^a + pQ)a[^a ~ po) 
^_ , B{po,a)da 



27r J2v 



2t: J2po a{\a + pQ)a{\a - poY 
where 

B{po,a) ^ [h{ia + p,) + h{-\a-pom\a-p)e-"'^- +h{-\a + p)e'^''-]. (26) 

The result is obtained in the following way. Wc assume that po > 0, so gets contribu- 
tions from 3 segments in the cr-axis, which are proportional to products of the ^'-matrix 
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elements, as follows: 



a e (-00, -2po) (-2po, 2po) (2po, +00) 



Tex S++S++ S++S__ (27) 



7^ oc — [-S-\ — S-\ — S-\ — S-\ — 5*- 



The arguments are p+ |cr in the first factor and p— |cr in the second factor. The S- matrix 
is given in Eq. (18), and the complex conjugation is taken on the real axis by means of Eq. 
(20), inverting the signs of the arguments. Thus the integral is given in terms of analytical 
functions with the general properties determined by the Schroedinger equation. 
Evidently, Lf is real and respects the reciprocity principle: 

Lt{x, xo) = L_t{xo, x). (28) 

Note that the result is translationally invariant; namely, for q ^ q — c the quantities r+ 
and a{K,) remain invariant, while ^ r_ — 2c and h{K) ^ 6(K)e^*''^. The second integral 
in TZ vanishes, if V{q) = V{—q), since in that case 6(k) = —h{—K) and B{p,a) = 0. The 
total transmission and reflection probabilities (for a given initial momentum p) are given 
by integration 

r+00 

T{p)= / r{p,r)dr=\a{p)\-\ (29) 



+00 



R{p) = / 7l{p,r)dr = \h{p) / a{p)\ 
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The integral representations in Eqs. (24-25) are manifestly causal. It is proven ||65| 



that the only singularities are due to zeroes of a(K), which are all in the lower complex 
half-plane for V{q) > 0. If r+ < 0, i.e. if q is ahead of the free propagation coordinate, 
the exponent for Im a > is decreasing, so the integration contour may be deformed to 
the upper half plane and the only contribution would be from 00 where a — 1, leading 
to T = 5(r+), as in free motion (8). For r+ > 0, the integral can be evaluated by means 
of residues in the lower half plane. Each zero of a{K) gives rise to a pair of poles in the 
integrand of (24) and (25). Thus T and TZ can be expressed as a sum over the ^-matrix 
singularities. The transmission propagator is given by, 

T(p, r+) = 5(r+) - ^(r+)S„Re {A^{p) exp[«2r+(p - /€„)]} , (30) 

where k„ are zeroes of a{n) and An{p) are determined by the corresponding residues, 9{r) 
is the step function (= for r < and = 1 for r > 0). Two examples are given below: 
the (5-potential barrier (Section 5) and the modified Poschl- Teller potential (Appendix 
C). In general, contributions from purely imaginary poles have the form of that for the 
5-barrier, Eq. (33). The functional dependence on the lag distance r_|_ is universal, though 
An and /€„ depend on the potential. As Im /t„ < 0, the second term in T{p, r^) is an 
exponentially decreasing and oscillating function of r+. It partly cancels contributions 
from the 5-function to the integrals for probabilities, so the effect of the barrier is to 
remove parts of the freely propagating wave packet. 
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5 Narrow potential barrier 

Let us consider the (5-potential barrier 

V\q) = Vo5{q), (31) 

for which the S'-matrix elements are given by 

a^K) = 1 - vo/2iK, b^n) = vo/2iK. (32) 

There is one pair of poles in the integrands in (24) and (25), and the integrals are calculated 
by the residue theorem, leading to 



T\p, r+) = 5(r+) - e{r+)2vo^l + {vo/Apye-''''+ cos(2pr+ + 7) (33) 



„;2 

'R.\p,r) = 9(r)^e-''°'- sin(2pr_ 
2p 



where 7 = arctan(fo/4p). The total transmission probability is 



T{p) = 1 



^0 



f + 4p2 

The 5-potential has a simple explicit solution also off the energy shell, 

^ ' ^ 27r(l + t;o/2v^) 

independently of k and ko. Thus one can evaluate the non-asymptotic time dependence 
ofEq. (14), 



+ 



ik' + \vl)iki + \v'o) 

1 (wiVitk^) w{^fM4 



(35) 



1.2 _ 1.2 \ 1.2 , 1„2 1.2 I 1,,2 



where W{z) is a function related to the probability integral p6 



W(z) = zw(z) = ze-'' (1 - erf(-iz)) x ^ f 1 + At + 0(z-^)] (36) 

V71" V 2z^ J 

In the large-time asymptotics, for t ^ max ^Vq^ , kQ^ , k"'^^ , jf vanishes as t"'^/^, more 
rapidly than in the general case since the transition matrix elements in Eq. (34) are 
insensitive to separation from the energy shell (cf. Appendix A). The result is 

As can be shown by means of the double Fourier transform in k and k^, the expres- 
sion in Eq. (35) agrees with the evolution kernel for the 5-potential in the coordinate 
representation, obtained previously[pO|. 
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6 Semi-classical approximation 



In the semi-classical approximation one has the eikonal formula, for the transmission 
amplitude 



a^n) = exp[iS{K)], S{k) = k 



+ CX) 



1 - VI - ''-'yii) 



dq. 



{31 



For > Vq = max[y(g)], 5* is real and \a\ = 1, so the transmission is complete. Note 
that S{—k) = —S{k), by virtue of the analytical continuation. Unlike the exact solution, 
the approximate function a(fi;) is not meromorphic, since S{k)/k has a branch point at 
= Vo. 

The resulting transmission propagator is given by 



T(p, r) 



TT Jo 



da cos 



(39) 



The total transmission probability is obtained by integrating T in r+. At large r+ the 
vicinity of a = is dominating in the integral, so one can use the series expansion in the 
exponent. For < Vq, 



S{p) + S{-p) = -2tl, I ^ I ^V{q)-p^dq, 



(40) 



where the integral is on the segment where V{q) > p^, / is real. That leads to the familiar 
semi-classical expression for the transmission probability, T = exp(— 2J). 

Unlike the integrated probability, the value of the propagator for a given r+ needs a 
more careful treatment. Expanding the exponent around a = and retaining terms up 
to the order of a^, one gets the Airy function for T{p, r+). While adequate for classically 
allowed regions 0], as well as for large r+, this naive semi-classical approximation can not 
be valid everywhere. In particular, it violates casuality in contrast to the results obtained 
from the exact quantum treatment. 

Fairly accurate results in some domains in the {pq, r_|_)-plane are obtained in the frame- 
work of the WKB approach if the integral in Eq. (39) is calculated by the stationary-phase 
method. If the critical points of the integrand appear in the region where the semi-classical 
approximation is valid, one can make use of Eq. (38). Positions of the critical points are 
given by the equation 



2r+ + S'i^a + p) + S'i^a - p) = 0, 

r+oo 



(41) 



V{q) 



dq. 



The latter has a classical meaning for > Vq; namely, S'{k) = —v{tv — tq), where 
v = 2ft: is the initial particle velocity, Ty and tq are the times of flight in presence of the 
potential barrier and in free space. Hence Eq. (39) can be interpreted in the following 
way. Treating a as a quantum fluctuation of the incident momentum, one can note that 
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r+ equals to the lag due to the barrier, averaged between two classical trajectories, their 
mean momentum being equal to p. For positive r+, the critical points of the exponent 
may take place at real u, and their positions depend on the magnitudes of r+ and p. 

At small r_|_ the region of large \a\ contributes substantially, and one use the high- 
energy asymptotics, where k — > oo, a(K) x 1, and 

5(«:)x^, w=r^V{q)dq. (42) 

Let us consider the domain where 

r+p > 1, p/w -C 1, r+p'^/w -C 1. (43) 

The inequalities mean that i) the domain is consistent with the uncertainty relation, ii) 
the process is a deep under-barrier tunneling, iii) the asymptotics of Eq. (40) can be used. 
The corresponding critical points are given by 

= 2w/r+ > (44) 

The contribution from the pair of critical points is 

r (p, r+) = - /°° cos (r+a + da x (^^^+) ^ ^os (2J2wr+ + ^] . (45) 
TT Jo \ a ) ^+v^ V 4/ 

This result is quite different from the standard semi-classical approximation, which is 
hardly surprising, since the domain given by inequalities (43) is not dominating in tun- 
neling as a whole, even though small r+ represent the most rapid signal transport. 

7 Summary and conclusion 

The phase-space propagator for a general local one-dimensional potential barrier has a 
universal form; it is given by a sum over singularities of the ^'-matrix. It conserves energy, 
is manifestly causal and shows certain features specific for quantum theory. 

At large times, the initial Wigner function is splited into three parts: the reflected bun- 
dle, the transmitted bundle and a transient group which can be neglected if the incident 
particle was prepared with a narrow momentum distribution. The resulting probabili- 
ties have a momentum dependence, which causes the known effects, like dispersion and 
forward attenuation. Asymptotically, after the particle leaves the interaction region, the 
time evolution of the Wigner functions is just the coordinate translation with constant 
classical velocities. As different parts of the Wigner function are translated with different 
velocities, the evolution of the phase-space distribution goes on, and the actual experi- 
ment results may depend on the detector position, in particular on its distance from the 
barrier region. 

The coordinate dependence of the propagators is universal for any local potential. 
The propagators are functions of the lag distance, which is the difference between the 
free motion and the current coordinate. At large lag distances, the probability has an 
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exponential decrease and is oscillating as a result of quantum interference effects. The 
decrease rate and the oscillation frequency are determined by positions of poles of the 
transition amplitudes in the complex energy plane. The coefficients and the phase shifts 
depend on details of the barrier shape and on the initial momentum. Of course, the 
phase-space propagator is not a (5-function. Each point in the initial Wigner function 
gives rise to disconnected domains in the final distribution. This is a purely quantal effect 
in the barrier penetration dynamics. 

The standard semi-classical approximation can be used for evaluation of the total 
transmission probability, obtained by integration of the wave packet outside the potential 
domain. It is hardly adequate, however, for description of the process time dependence, 
since the analytical properties representing the causality in the energy representation are 
not maintained properly. 

In contrast to classical theory, individual phase-space trajectories cannot be traced in 
quantum theory. Both the the initial distribution and the detector acceptance have final 
supports, and the observed probabilities are obtained by integration. It is evident, how- 
ever, that the principle of causality is not violated by quantum theory, and no information 
can be transported faster because of the barrier. The problem of the signal transport and 
the time delay is a subject of a forthcoming paper |]67|. 
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Appendix A 

OfF-energy-shell corrections in the large-i asymptotics 

The transition operator of Eq. (11) has the following general property 



which leads, in particular to the ^-matrix unitarity. Hence one has for the discontinuity 
of the matrix element in Eq. (14) 



(46) 




TT 



(47) 



D{K;k,ko) = {k\fs\uK){ko\fs\uK), 



u=±l 



where n = e^l"^ . As a function of k, D is regular. Besides, it has the following general 
property of positive definiteness. 





for any integrable complex function /(/c). In particular, -D(k; k,k) > for all k. 
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The "half-on-shell" transition amplitudes present here are expressed in terms of solu- 
tions of the Schrodinger equation (19) y±{x), satisfying the complimentary asymptotical 
conditions, as given in Ref . [|65[] , 

{k\TAvK) = / rfxe-'(^-"")"|(/t2-A;2)(r/U_-r/+r7M (49) 

+\/(x) [((2 + v)k - k)y+{x)e-'^'' + ((2 - u)k + k)y.{x)e+'^''] } , 

where t]±{x) = y±{x) —e^^'^^ for x — * ±oo. It is evident from the integral representa- 
tion that the transition amplitudes are regular functions of k, and they are proportional 
to K at /c^ = K^. (Note that na is finite at k = 0.) 

Finally, the off-energy- shell correction is given by the following equality (we assume, 
for simplicity that there are no bound states) 

.Mk,k„) = exp|it(i' + kl)] p.v.^~e-'^-5^<|AM^ ^ O(rV^). (50) 
The large-t asymptotics of the integral is determined by the behavior of D at k ^ 0. 



Appendix B 



Gaussian distributions 

Let us consider, for example, the Gaussian initial phase-space distribution po{q,p) and a 
similar detector acceptance function Cil^p)- We shall assume, for simplicity, that they 
describe pure states which have the minimal uncertainties, allowed by the Heisenberg 
principle. The functions are 

^VAo-Ao(p-Po)'] , 



exp[-(g - (5o)7Ao 



exp[-(g - Qd) /Xd 



(51) 



where {Qo, Pq) and {Qd, Pd) are the most probable coordinate and momentum values 
for the initial wave packet and the detector, respectively; Aq and A^ are the coordinate 
dispersions. If the wave packet was prepared with a fairly definite initial momentum Pq, 
then Ao must be relatively large, while A^ is small for a detector localized in space. 

The probability to detect the particle at a time t is given by Eq. (5). For pure states 
the probability is expressed in terms of the transition amplitude. 



Mt) = -(AoAd)'/'exp 



AnA 



An + A, 



\A + B\ 



where 



A 



dn 



B 



exp 



-i(A„ 



Xd){K- P+y -Ik{Q_ + Kt) 



(52) 



(53) 



oo a K 



6(k) exp 



P± = (AoPo±AdPd)/(Ao 



P_)2 
-Xd), 



Q± = Qo^Qd 



13 



Apparently, if the ^-dependence of \a\ is not strong, the region of k ~ P+ is the most 
substantial in the integral for the transmission amplitude A{t). On the other hand, the 
amplitude vanishes in the large-t asymptotics except for times t given by the stationary 
phase condition, 

vt = Qd - Qo - dip/ dn, (54) 

where w = 2k is the velocity in the substantial region, and <^(k) = arg(a). Similar 
arguments can be applied to B{t). (A more detailed discussion is given in Ref.|]67|.) 

The probability in (52) is a sum of three terms, = Wt + Wr + '2ws, corresponding 
to transmission, reflection and the transient interference, as in Eq. (23) for Lj. At large 
times t and the states localized outside the potential region (large Qo and Qj), one of 
the amplitudes is exponentially smaller than the other one. As soon as w"^ = WfWr, 
the interference term can be neglected. In general one deals with mixed states, where 

< WtWr, so the interference is always negligible. 



Appendix C 

Modified Poschl- Teller potential barrier 

The modified Poschl- Teller potential barrier, 

V'^^iq) = (55) 
cosh [q/s) 

where Vq and s are constants, is an example for a local potential barrier which enables an 
explicit solution and a test of the various approximations. The known solution is 



SKr(| + u + isfi;)r(| — uj — isn) 



Smh TTSK 



where u = — tigS^; it is real if the barrier is narrow. 

The singularities of the transmission amplitude are simple poles at sk = —i{n + j±Lj), 
n = 0, 1, ■ ■ ■ . The phase-space propagator is given by an infinite sum of the residues. For 
large t the result is 

L^'^{q,P',Qo,Po) >^S{p-po)T^'^ir+,pQ) +5{p + po)n^^{r_,po), (57) 
rPT(r+, po) = S{r+) + 9{r+) \Tt{y, uj) + Tt{~v, uj) + Tt{y, -uj) + Tt{-v, -uj)\ , 
7^^'^(r_,po) = Q{rJ) \Tr{v-,uj) + :F^(-z/, cj) + :Fr(z/, -cj) Ar T r{-v , -uj)\ 
+^(-r_) \T,{v, uj) + T,{-v, -uj)\ 

where v = 2pqs and, 

2r{2uj)r{iu)r{iu + 2uj) 
'^""'"^^ ~ sTi-l + uj)ra + a;)r(-i + + u;)ra + ZU + U) 
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exp[(iz/ + 2uo- 1)—] 4F3 



— uo,\ — uj,\ — iv — uo,\ — iv — uj 
1 — 2uj, 1 — ivA — iv — 2uo 



exp(-2-5^) 
s 



s[r(i + ^)]2r(i - ^)r(-i + a;)r(-i + zz/ + 



exp[(iz/ + 2ct; — 1) — ] 4F3 



I - oj, i - a;, I - zz/ - ^, i - zz/ - a; _ ,2!:: 
1 - 2tu, 1 - zz/, 1 - iz/ - 2a; ' s 



r(f + W + cu)r(f + uj)Y{-l - iu) 
Y{\ + uj)Y{\-uj)Y{l + iv) 



exp[(iz/ + 2)— ] 4F3 



I + CJ, I - tU, I + ZZ/ + Co", I + zz/ - w 

2, 1 + iz/, 2 + iz/ 



,2r_ 



exp( 



(5^ 



The 4F3 are the generahzed hypergeometric functions |]68 

ill ^2-1 i'ii 



4-^3 



-^1; -^2, A3 



c 



C (6)n(6)n(e3)n(e4) 
^0^! (Al)n(A2)„(A3)„ 



(59) 



and {i)n = r(^ + n)/r(^). The series are convergent, as long as \(\ < 1. 

For a wide barrier, uj is imaginary, and the series approach their radius of convergence 
as s ^ 00. Thus the integral representation is more practical than the series. For small 
momenta, k <ti Vq, one gets: 



a (k) ^ exp[—TTs{vo — k) + iS{K)], 
= SK,m[— — 1) + sfoln ■ 



(60) 



Vq- K 



which is the semi-classical approximation (the integral in Eq. (38) is calculated). Now T 
is given by the Airy function. 



PT 



-2.s(.o-Po)(3^3)-i/3^^ [(3^3)-'/'(r+ - ^1) 



(61) 



V^i = -sin -| 



- 1 



V^3 



The appearance of the Airy function is typical for semi-classical Wigner's functions. As 
s ^ 00, the Airy function approaches the 5-function, corresponding to the classical prop- 
agation, 



FT 



-27rs(t;-po)^ 



(62) 



The additional term looks like an advance in time, violating causality, yet in fact it is an 
artifact of the naive semi-classical approximation. 

In the opposite limit of a narrow barrier, i.e. small s, the series converges rapidly. If 
we set s — > 0, VqS =const, the 5-function discused in Section 5 is reproduced. One has no 
problems with the causality in that approximation. 
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